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Strongly coupled quantum heat machines 


David Gelbwaser-Klimovsky and Alan Aspuru-Guzik 
Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138 

Quantum heat machines (QHMs) models generally assume a weak coupling to the baths. This supposition 
is grounded in the separability principle between systems and allows the derivation of the evolution equation 
for this case. In the weak coupling regime, the machine’s output is limited by the coupling strength, restricting 
their application. Seeking to overcome this limitation, we here analyze QHMs in the virtually unexplored strong 
coupling regime, where separability, as well as other standard thermodynamic assumptions, may no longer hold. 
We show that strongly coupled QHMs may be as efficient as their weakly coupled counterparts. In addition, we 
find a novel turnover behavior where their output saturates and disappears in the limit of ultra-strong coupling. 


One of the basic tenets of standard thermodynamics is the 
principle of separability, which allows to clearly define and 
distinguish systems that interact with each other. When the 
surface to volume ratio is small, surface effects are negligible, 
and thermodynamic variables only depend on the volume and 
not on the shape. This argument implicitly assumes a weak 
coupling, restricting the interaction space to a small interface 
between the systems 

The assumption of weak coupling was essential for the de¬ 
velopment of open quantum system theory a, in particu¬ 
lar for the development of the Kossakowski-Lindblad master 
equation aa, that describes the evolution of a system inter¬ 
acting with a thermal bath. Quantum heat machines (QHMs) 
models EHIl use this framework to describe the evolution 
of the “working fluid” under the influence of the hot and 
cold baths. Progress in this field has been recently reviewed 
insiiiii. QHMs may operate either as engines, by extracting 
work power, or as refrigerators, by investing work power and 
cooling the cold bath. In both cases, quantum resources have 
been proposed lfT7l420l in order to boost their output and ef¬ 
ficiency. Nevertheless, these models assume a weak coupling 
to the baths, resulting in limited QHMs outputs and conse¬ 
quently restricting their applications. 

The potential technological implications of high-output 
QHMs, such as faster and more powerfull laser cooling 
iiziiiaa, call for a prompt way to overcome the limitation 
set by the weak coupling assumption. However, the strong 
coupling limit has been virtually left unexplored due to the 
lack of theoretical tools to describe the “working fluid” evo¬ 
lution. One of the few exceptions ll23l considers the case of 
Hamiltonian quench, which involves the switching “on” and 
“off” of the system-bath interaction Hamiltonian, introduc¬ 
ing an energy and efficiency cost that reduces the machine 
efficiency below the Carnot bound. In this letter we take a 
different approach by putting forward a strongly coupled con¬ 
tinuous QHM model (see Fig. [T]i, that does not require the 
coupling to and uncoupling from the baths, which may not be 
possible at nanoscale, where the system is totally embedded 
in thermal baths. We investigate its output and efficiency in 
order to determine its performance limits and which thermo¬ 
dynamic principles, e.g., Carnot bound, still hold at the strong 
coupling regime. Addressing these issues becomes more rel¬ 
evant in the light of the large progress achieved in the field of 
strongly coupled superconductors Il24l427l . which makes the 
realizations of strongly coupled QHMs potentially tractable 


in the near future. 
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Figure 1: Model of continuous quantum heat machine where the cold 
bath strongly interacts with the working fluid, while the hot bath is 
weakly coupled. 

Model and analysis We employ a model for a continuous 
QHM similar to the one we studied previously in the weak 
coupling limit ll28l . This system can operate either as an en¬ 
gine or a refrigerator depending on the spectrum of the reser¬ 
voirs and the engine’s driving frequency. This model is com¬ 
prised by a driven two-level quantum system, that represents 
the working fluid, permanently coupled to the heat baths (hot 
and cold). The evolution of this model is governed by the 
Hamiltonian 

n = + 2(a+e-“'‘ + a_e“'*)+ 

o-z O X ^c{gc,ka\ + 5c,fcOfe) + 

k 

CTx O X ^H{gH,khl f-gH,kbk) 

k 

+ '^ujc,kalak -\-^ujH,kblhkj ( 1 ) 

k k 

where ^c{h) is the strength parameter of the cold (hot) bath, 
Pi fe is a dimensionless parameter that defines the relative cou¬ 
pling strength of the TLS to the mode k of the i-bath, ay are 
the standard Pauli matrices, and a\, ak {b\,bk) are the creation 
and annihilation operator of the cold (hot) bath mode k. The 
election of the hot and cold bath is somehow arbitrary and a 
similar analysis could be performed if they are interchanged. 
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The coupling is consider weak if 'jTcor 'C 1, where 7 is 
the decay rate and is equivalent to the resonant coupling spec¬ 
trum (7 = G{uJo)) and Tcor is the bath correlation time Il29l . 
Is it possible to extract work or to cool down in the strong 
coupling regime? To elucidate this question, we consider that 
both couplings are strong. While the reduced dynamics is an¬ 
alytically solvable in the weak regime, in this case the per¬ 
turbation expansion on the coupling strength contains infinite 
no-neglectable terms lf30ll . 

Nevertheless for H, this obstacle may be overcome by solv¬ 
ing the problem in a more appropriate basis, where the sys¬ 
tem is effectively weakly coupled to the two baths. This is 
achieved by using the polaron transformation EDEl, e®, 
where S = az and ak = 

The transformed Hamiltonian, "H = is 


77 Wo 

H 






^ 0 {A+ -A)+ ® (^_ - A)) + 

(cr+ (g) A+ -f tT_ (g) y4_ ) (g) E ^H{gH,kb\ + g*H,k^k) + 

k 

ujc,ka\ak + ujH,kb\bk, ( 2 ) 


where A± = XlkD{±2ak), D{ak) = 

is the displacement operator, A = (^±) = 


coth(- 




e ■'-““II “-fc II ' ^ qqjg terms 

on the Hamiltonian proportional to the identity have been 
neglected. 

In the transformed Hamiltonian, the coupling operators are 
different. A+ — A and A- — A, instead of a \, ak and ex¬ 
tra terms are added to the hot bath coupling (see Eq. © and 
Suppl. A). As we show below, the new couplings may be ef¬ 
fectively weak even for high values of the original coupling 
strengths, Cc(rr)- Therefore, the assumptions derived from 
the weak coupling are correct (e.g. the transformed baths re¬ 
main at thermal equilibrium) and the master equation may be 
derived using standard techniques l|4] [38l also for values of 
^C(H) that break the weak coupling assumption in the origi¬ 
nal basis ll32H37]l . 

The transformed cold bath, imw interacts with the TLS 
through two different operators, Fi {t) = ^ (A_ (f) — A) and 

F 2 = A_{t) (g) Y.k {gH,kbi{t) + gH,k^k{t)Y The corre¬ 
lation function of the first is 


(Fi(f)tFi(0))=(^ 




Afc(t) 


- 1 ), ( 3 ) 


where 
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Figure 2: (Color online) Effects of the coupling strength. Main panel: 
Coupling spectrum in the original basis Gi(a;o) (dotted line) and in 
the transformed basis Gi(a;o) (continuous line) as a function of the 
coupling strength ~ Crr- The weak coupling assumption holds 
only for coupling spectrum below the dashed line. Inset: Power for 
a QHM as a function of the coupling strength, ~ A turnover 
is observed and the power decays for ultra-strong coupling. At this 
limit, the lack of system-bath separability prevents work extraction. 


ecY^kkit) = {Fi{t)Fm) = 

k 

E^cliffc.fef ^cos(a;fcf)coth(^^^^)-fsin(a;fcf)^ (4) 

is the time correlation of the original coupling operator, 
Pi{t) = EfcCc (sc.fcaiWand /3c is the 
equilibrium temperature of the transformed cold bath. The 
coupling spectra that govern the evolution are derived from 
the correlations of the transformed operators, Gi{uj) = 

zG 1 , 2 . 

In Fig. the dependence on the coupling strength, 
^c^ of the coupling spectrum in the original = 

dotted line^ and transformed ba¬ 
sis ^Gi(a;), continuous line^ are compared. While both cou¬ 
pling spectra are proportional to the square of the coupling for 
small coupling strengths, in other regimes their behavior di¬ 
verge. The validity of the “weak” coupling assumption for the 
spectrum Gi(a;) has been broadly shown Il32l437l l^ . In a 
similar manner, the operator A±{t), will constraint G^ioj) to 
the weak coupling regime as long as ^c- 

Gi(a;) keeps the standard KMS condition Gi(—w) = 
EiEIl. It includes modes harmonics, i.e., 
Gi(a; > LOcutoff) 0 as long as w is a linear combination of 
bath modes harmonics. This propriety lets the use of highly 
detuned baths in strongly coupled QHMs, unlike for weakly 
coupled QHMs that require resonant baths (or at least resonant 










3 


with linear combinations of uiq and w;, the TLS and driving 
frequency, respectively HI). 

The polaron transformation allows the derivation of the 
QHM evolution for a wide range of values of the coupling 
strengths. Nevertheless, this simplihcation entails other com¬ 
plications, as the loss of separability. In the transformed ba¬ 
sis, the second correlation is far from standard. It involves 
exchange of excitation with both baths (the operators and 
A± for the hot and cold bath respectively).The lack of sepa¬ 
rability breaks the standard KMS condition, casting doubt on 
the validity of other thermodynamic principles, as the Carnot 
bound. 

The answer to this question is obtained from the theory of 
non-equilibrium thermodynamics, which introduces the fre¬ 
quency dependent “local” temperatures, /3(w). They are anal¬ 
ogous to the non-equilibrium position-dependent local tem¬ 
peratures ll42l . In the non-equilibrium framework, the KMS 
condition is generalized (see Suppl. B): 


G2(-cc) = e-«‘")‘"G2(cc), 

/3(u;) =/3cA(a;)-f/3ff (1-A(a;)), (5) 


where A(a;) measures the relative contribution of the trans¬ 
formed cold bath to G 2 (uj) and may take any positive or neg¬ 
ative value. Therefore j3(u!) is not restricted to the range 
[PhtPc]- It depends on both baths coupling strength distri¬ 
bution and modes, making /3(a;) frequency dependent, blur¬ 
ring its physical interpretation. As we show later, it allows to 
establish clear thermodynamic bounds to the efficiency of the 
QHMs and to relate them to the Carnot bound. The precise 
value taken by /3(a;), depends on how the exchange energy w 
is divided between the hot and the cold baths. 

In a similar manner as Gi (w), G 2 (w) not only includes har¬ 
monics of the cold bath, but combinations of them with modes 
of the hot bath. Therefore, also the hot bath may be highly 
detuned from the TLS frequency (and from any linear combi¬ 
nation with the driving frequency). 

In the transformed basis, we use the standard weak coupling 
master equation based on the general Floquet theory of open 
systems ||38l. We just stress the main steps, but the detailed 
derivation may be found in ll28l . The reduced evolution of 
the TLS density matrix, p, is given by a linear combination of 
Lindblad generators obtained from the Fourier components in 
the interaction picture of the working fluid coupling operators 

and = 

SqGZ Sw 'S'2,g(w) e In the interaction picture. 


dp 


^=cp. c=y.a+j:^ 


2 

qoj 


q,uj q,uj 

m G\uj + quJi)^^ 
L-qpp- 2 "" 


( S,^q{^)p,Slg{uj) + S,^q{uj),pSl^{uj) 


( 6 ) 


The TLS density matrix evolves until it reaches a steady 


state (limit cycle), Cp = 0. At this point any transient ef¬ 
fect averages out and one may calculate the steady state work 
power and heat flows. 


Ji = ^ sgn(a;) (w -f quJi) Tr [L\^p\ , P = -Ji- J 2 , 

quj 

(7) 


where sgn(a;) = 1 for oj > 0 and sgnjw) = — 1 for w < 0. 

In particular we are interested in the ultra-strong coupling 
regime ~ ^ ^ 0 ) to And out if QHMs may have an 

ultra-high output. Nevertheless at this limit, and assuming a 
weak driving with positive detuning ((5 = ojQ—ijJi ^ > 0), 

the work power dependence on the coupling strength goes as 
(see Suppl. C) 


P oc oc 

sc 




oth( ) 


SC 


( 8 ) 


The conditions for work extraction {P < 0), 


^ < 1 _ ^(^ 0 ) 

wo Pc ’ 

is derived from Eq. The heat currents to the baths are: 


Ji cx < 0, 

sc 

J 2 a > 0. (9) 

sc 

For ultra-strong coupled baths, the work power will decay 
with the coupling strength as shown in Eq. The exact 
counterpart of Eq. ([^ for any value of ^c ~ is plotted on 
figure l^lnset. Opposite to what may be expected from pre¬ 
vious results in the weak coupling regime, work power does 
not increase indefinitely with the coupling strength. Not only 
it saturates, but at some point, decays and vanishes. At the 
ultra-strong limit, the system and the baths are no longer in¬ 
dependent, preventing work extraction which requires some 
degree of separability. 

The determination of the engine efficiency, as well as the 
cooling power (in the refrigerator operation mode), is a more 
subtle issue. A naive guess would be to consider J 2 , which is 
positive, as the incoming heat flow from the hot bath and to 
define the efficiency as 


-P Wi P{oJo) 

t? = ^ = — < 1-3—> 

J 2 Wo Pc 


( 10 ) 


which can take any value, even above Carnot limit. Never¬ 
theless, the lack of separability complicates the determination 
of how much energy is exchanged with each bath through the 
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coupling spectrum G 2 (w). Only a fraction, 1 — A(a;), of J 2 
is originated in the hot bath. Therefore the correct efficiency 
expression is 

— _ tU; Ph _ 

^ J 2 (1 — A (wo)) Wo (1 ~ A(wo)) Pc 

( 11 ) 

From Eq. ([TT) we conclude that the Carnot bound may 
be reached, but not surpassed, by the appropriate choice of 
driving frequency. Therefore, strongly coupled machines are 
as efficient as their weakly coupled counterpart. 

In a similar way, one can calculate the cooling power for 
the refrigeration operation. This sets the opposite condition 
on the frequencies; ^ > 1 — ^ making Ji > 0, which 

can erroneously be confused with the cooling power. The lack 
of separability between both baths (the dependence of on 
both baths operators) mixes the heat flows between both baths 
and part of J 2 is heat flowing to the cold bath. The correct 
expression for the cooling power is Jq = Ji — A(wo) J 2 and 
is limited by the Carnot bound for refrigerators. The cooling 
power has a similar dependence on the coupling strength as 
the work power and also decays and vanishes for ultra-strong 
coupling. 

An ideal platform to test our results are superconducting 
quantum circuits, where the almost unexplored strong cou¬ 
pling regime has been recently achieved Il43ll44ll . showing as¬ 
tounding ^i/ujo ratios of around 0.12. Moreover, recent theo¬ 
retical studies have shown that a (Ta;—coupling between quan¬ 
tum microwaves and artificial Josephson-based atoms can be 
pushed up to ^i/ujQ ^ 2 ll45l . which is well beyond the crit¬ 
ical point (^i/ujo = 1) at which the power efficiency is max¬ 
imum. Our proposal consists of a periodically driven super¬ 
conducting flux qubit with tunable gap ll46l . where the main 
loop is coupled to the hot bath (cx coupling), and the a—loop 
is coupled to the cold bath (a^ coupling). In order to bring the 
Cz coupling to the strong regime, we galvanically couple the 
a—loop to the open transmission line that plays the role of the 
cold bath. 

Conclusions The possibility of work extraction and cool¬ 
ing in the strong coupling regime was shown. Even though 
some thermodynamic principles, as the standard KMS con¬ 
dition, do not longer hold at this regime due to the lack of 
separability between the baths, the operation of the QHMs 
may be described in a non-equilibrium framework. This is ad¬ 


vantageous because it shows that important principles, as the 
Carnot bound, still hold in the strong coupling regime. The 
introduction of frequency-local temperatures, that account for 
the different baths contributions to the heat flows, are useful 
to determine how the heat flows are divided between baths 
and to correctly calculate the QHMs efficiency. As we have 
shown, continuous strongly coupled QHMs, as their weakly 
coupled counterparts, avoid the efficiency reduction due to the 
coupling turning on and off and keep the Carnot bound which 
can be reached under the appropriate driving frequency. 

The appearance of the “non-equilibrium” temperatures is 
related to the loss of separability. Even though both baths, 
in the transformed basis, are in equilibrium, the heat flows 
mixes the contribution of both of them, causing an effective 
deviation from equilibrium. 

There are similarities between weakly and strongly coupled 
QHMs, but the differences should not be overlooked. While 
weakly coupled QHMs require baths with modes resonant to 
the TLS (or linear combinations with the driving frequency), 
strongly coupled QHMs operate also for highly detuned baths, 
because harmonics of the strongly coupled bath modes also 
contribute to both coupling spectra. An important feature of 
strongly coupled QHMs is that, differently to their weakly 
coupled counterpart, where the outputs are proportional to the 
square of the coupling strength, work and cooling power sat¬ 
urate at some point and for ultra-strongly coupled machines 
they fall down as the coupling strength increases. This is a 
consequence of the lost of separability as the coupling strength 
increases, and shows that QHMs require some degree of sep¬ 
arability to operate. 

In order to optimize QHMs output the “right” coupling 
strength is needed, resembling the quantum Goldilocks effect 
BtII found in photosynthetic systems. The latter should be fur¬ 
ther investigated to determine if evolution fine-tuned the cou¬ 
pling strength to the baths in order to maximize their chem¬ 
ical power output. Alternatively, the turnover behavior may 
be corroborated experimentally using superconducting qubits 

mM- 
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Supplementary Information 

A. System-Bath coupling 

In the original basis the system-bath coupling operators are; 


® ^ ic{.gc,kal + gcM^k), 

k 

crx<S>Y^ ^H{gH,kb\ + gH,kbk), 

k 

where ai are Pauli matrix and operate on the system. a\. and ( 6 j. and b/,) are the cold (hot) bath operators. 
In the transformed basis, the system-bath coupling operators are 

^ (g) {A+ -A) + e“'V_ 0 {A_ - A )), 

(cr+ (g) A+ -I- cr_ (g) A_) (g) E (HigH,kbl-\- gH,kh), 

k 

where 


(SI) 


(S2) 


A±-A = 


oth(- 


(S3) 


B. Generalized KMS condition 

As mentioned in the main text, the coupling spectrum G 2 (w) contains contributions from both baths. The frequency sum of 
the contributing hot and cold bath modes should match the spectrum frequency, uj = ujH,i + There are many combinations 
of modes that match the spectrum frequency, therefore 


G2{uj) — G2{uJH,i + ^C,j)- (S4) 

Due to the non-linearity of the cold bath coupling operators in the transformed basis, its mode harmonics also contribute to 
the sum on Eq. ( |S4| l. The G 2 {oJH,i + ^C,j) physical meaning is an energy exchange, where an excitation uj of the system is 
interchanged with the LOH,i, and tucj modes of the hot and cold baths, respectively. They keep a modified KMS condition 

G2i-iOH,^ - OJCJ) = + wcj). (S5) 

Combining all the terms, the effective frequency-local temperature may be dehned as 


^-P{u))uj 


G2{-UJ) 

G2{uj) 


E If TfT 

i,3 


G2{^H,I + ^C,j) 
G2{i^H,i + ^C,j) 


where ^ is the relative weight of the G 2 {ojh,i + ^C,j) component. 


(S 6 ) 


C. Heat currents and power 


For a weak driving and positive detuning, the heat currents and power are (see ESlH 
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J2 = —Wo 


Gi((5) + G2(wo) 
Gi(<5)G2(wo) 


— e 


-/3((^o)‘^o'| 


Gi(,5) + G 2 (wo) 


(e 


-Pc& _ g-/3(‘^o)‘^o'| 


P = UJi 


jGi(J)G^(wo) 

Gi((5) + G 2 (wo) 


,-PcS _ -P(u)o)i^o 


(57) 

(58) 

(59) 


where G,(w) = e*“"(F,(t)tF,(0))dt. 

We assume that the main contribution to the coupling spectrum comes from few modes. For the sake of simplicity we present 
the calculation assuming that this contribution is due to one mode. Then, 




Afc(t) 


2 


- 1 ) ~ ^ j (e '“ 0 - 1 ) 


nA 

~Y 


2/00 


g 2 J 2 „ ( ] Cos{2nuJk,t) + 2iJ2n+i ( ) Sin{{2n + l)wfc„f) 


\n—Q 




n =0 


^0 / 


^ 2/„ 4 




coth( ^^^^° ) \ Cos(nujkot) 


(SIO) 


\n—0 


^ko 


) 


where Bessel and modified Bessel functions have been used to expand the exponential. Using the Fourier transformation and 
taking the asymptotic limits of the Bessel and modified Bessel functions: 


-4eci:,\\"-^\\" cothi^) 

Gi(<5)(x --2 -. (Sll) 

sc 

G 2 (wo) has a similar dependence. Therefore, for —> 00 , P (x - 72 -• 



